Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
伊藤 孝; 門野 良典*
Journal of the Physical Society of Japan, 93(4), p.044602_1 - 044602_7, 2024/04
We propose a model to describe the fluctuations in the internal magnetic field due to ion dynamics observed in muon spin relaxation (
SR) by an Edwards-Anderson-type autocorrelation function that separates the quasi-static and dynamic components of the correlation by the parameter 
(where 
). Our Monte Carlo simulations for this model showed that the time evolution of muon spin polarization deviates significantly from the Kubo-Toyabe function. To further validate the model, the results of simulations were compared with the SR spectra observed in a hybrid organic-inorganic perovskite FAPbI
[with FA referring to HC(NH
)], where local field fluctuations associated with the rotational motion of FA molecules and quasi-static fields from the PbI lattice are presumed to coexist. The least-squares curve fitting showed reasonable agreement with the model with 
, and the fluctuation frequency of the dynamical component was obtained. This result opens the door to the possibility of experimentally distinguishing fluctuations due to the dynamics of ions around muons from those due to the self-diffusion of muons. On the other hand, it suggests the need to carefully consider the spin relaxation function when applying SR to the issue of ion dynamics.
Riyana, E. S.; 奥村 啓介; 坂本 雅洋; 松村 太伊知; 寺島 顕一; 神野 郁夫
Journal of Nuclear Science and Technology, 8 Pages, 2024/00
被引用回数：0 パーセンタイル：0.18(Nuclear Science & Technology)We investigated the possibility of estimating the effective neutron multiplication factor (
) of the fuel debris inside the canister and primary containment vessel (PCV) of Unit 2 of the Fukushima Daiichi Nuclear Power Station (1F) using remote gas-radioactivity measurement via simulation-based calculations. Our results demonstrate an almost linear correlation between and the 
Kr-to-
Xe activity ratio with respect to various fuel debris compositions. This correlation is maintained regardless of geometries such as the fuel debris canister and the PCV.
Tuya, D.; 長家 康展
Journal of Nuclear Engineering (Internet), 4(4), p.691 - 710, 2023/11
モンテカルロ中性子輸送計算手法は、固有値や積分中性子束などの様々な量を正確に評価するために用いられる。しかし、ある分布量を求める場合、モンテカルロ法では連続的な分布が得られないのが一般的である。近年、モンテカルロ法で連続分布を得るために、関数展開法やカーネル密度推定法が開発されている。本論文では、モンテカルロ法によって得られた訓練データと人工ニューラルネットワーク(ANN)モデルを用いたある物理量に対する連続分布の推定手法を提案する。概念実証として、2つの核分裂体系における反復核分裂確率(IFP)の連続分布をANNモデルにより推定した。ANNモデルによるIFP分布を、元のデータセット及びPARTISNコードで得られた随伴角中性子束と比較した。比較の結果、一致や不一致の程度はさまざまであったが、ANNモデルはIFP分布の一般的な傾向を学習することを確認した。
涌井 隆; 高岸 洋一*; 二川 正敏
Materials, 16(17), p.5830_1 - 5830_16, 2023/09
被引用回数：0 パーセンタイル：0(Chemistry, Physical)水銀ターゲット容器は、陽子ビームの入射に伴い、キャビテーション損傷を受けるため、キャビテーションバブルの位置や衝撃圧力分布の不確実さを考慮して、モンテカルロ・シミュレーションを用いたキャビテーション損傷を予測する手法を提案した。本手法では、個々のキャビテーション気泡の崩壊に起因する衝撃圧力の分布はガウス分布とし、衝撃圧力の最大値の確率密度分布は3種類の分布: デルタ関数、ガウス分布、ワイブル分布と仮定した。衝撃圧力の分布を記述する方程式の2つのパラメータについて、実験から得られたキャビテーション損傷の分布とシミュレーションから得られた累積塑性ひずみの分布を比較し、ベイズ最適化を使用して推定することができた。また、ワイブル分布を用いて得られた結果が、他の結果に比べて、実際のキャビテーションエロージョン現象を再現することが分かった。
久保 光太郎; 田中 洋一; 白田 勇人*; 荒毛 大輔*; 内山 智曜*; 村松 健
Mechanical Engineering Journal (Internet), 10(4), p.23-00051_1 - 23-00051_17, 2023/08
福島第一原子力発電所事故後、原子力発電所の外部事象に対する確率論的リスク評価(PRA)の重要性が再認識された。地震PRAでは、システム、コンポーネント、構造(SSC)の相関故障の扱いが非常に重要である。SSCが連関して故障した場合に、事故緩和システムの冗長性に悪影響を与える恐れがある。原子力機構では、SSCの相関故障を詳細かつ現実的に取り扱うために、Direct Quantification of Fault Tree by Monte Carlo Simulation(DQFM)法と名付けられたフォールトツリー定量化手法を開発している。この手法では、地震応答や耐力に関する相関を有する不確実さをモンテカルロ・サンプリングで考慮することにより、頂上事象の発生確率を定量化することができる。DQFMの有用性は既に実証されているが、その計算効率を向上させることで、リスク分析者は複数の分析を行うことができるようになる。そこで本研究では、簡易的な地震PRAのDQFM計算に準モンテカルロ法および重要度サンプリングを適用し、少ないサンプリングで分岐確率を定量化する効果を検証した。準モンテカルロサンプリングを適用することで、中・高地震動レベルでの結果の収束性がモンテカルロサンプリングに比べて一桁加速されることを示すと共に、重点サンプリングの適用により、低い地盤変動レベルにおいて故障確率が低いSSCに対して統計的に有意な結果を取得できることを示した。これらのサンプリング手法の応用による計算の高速化はリスク情報を活用した意志決定におけるDQFMの実用化の可能性を高めるものである。
Li, C.-Y.; Wang, K.*; 内堀 昭寛; 岡野 靖; Pellegrini, M.*; Erkan, N.*; 高田 孝*; 岡本 孝司*
Applied Sciences (Internet), 13(13), p.7705_1 - 7705_29, 2023/07
被引用回数：1 パーセンタイル：0(Chemistry, Multidisciplinary)For a sodium-cooled fast reactor, the capability for stable cooling and avoiding re-criticality on the debris bed is essential for achieving in-vessel retention when severe accidents occur. However, an unexploited uncertainty still existed regarding the compound effect of the heterogeneous configuration and dynamic particle redistribution for the debris bed's criticality and cooling safety assessment. Therefore, this research aims to develop a numerical tool for investigating the effects of the different transformations of the heterogeneous configurations on the debris bed's criticality/cooling assessment. Based on the newly proposed methodology in this research, via integrating the Discrete Element Method (DEM) with Computational Fluid Dynamics (CFD) and Monte-Carlo-based Neutronics (MCN), the coupled CFD-DEM-MCN solver was constructed with the originally created interface to integrate two existing codes. The effects of the different bed configurations' transformations on the bed safety assessments were also quantitively confirmed, indicating that the effect of the particle-centralized fissile material had the dominant negative effect on the safety margin of avoiding re-criticality and particle re-melting accidents and had a more evident impact than the net bed-centralized effect. This coupled solver can serve to further assess the debris bed's safety via a multi-physics simulation approach, leading to safer SFR design concepts.
涌井 隆; 高岸 洋一*; 二川 正敏; 田邉 誠*
実験力学, 23(2), p.168 - 174, 2023/06
核破砕中性子源の水銀ターゲット容器において、陽子ビームが水銀に入射することにより、容器内面にキャビテーション損傷が生じる。気泡核の位置や衝撃圧力分布のばらつきを考慮して、モンテカルロ・シミュレーションを用いたキャビテーション損傷の予測方法を提案した。実験より得られたキャビテーション損傷の分布とシミュレーションより得られた積算ひずみの分布を比較し、ベイズ最適化を用いた逆解析により、衝撃圧力の分布を評価した。ガウス分布を仮定した最大衝撃圧力の平均値及び広がりはそれぞれ3.1GPa及び1.2mであった。シミュレーション結果は実験結果を再現しており、本評価手法が有効であると言える。
Li, C.-Y.; 渡部 晃*; 内堀 昭寛; 岡野 靖
Proceedings of 30th International Conference on Nuclear Engineering (ICONE30) (Internet), 10 Pages, 2023/05
Quantitative assessment of the effect of accident management on the various external hazards is essential in the nuclear safety analysis. This study aims to establish the dynamic probabilistic risk assessment methodology for sodium-cooled fast reactors that can consider the transient plant status under continuous external hazards with corresponding countermeasures operating stochastically. Specifically, the Continuous Markov chain Monte Carlo (CMMC) and Deterministic and Stochastic Petri Nets (DSPN) methods are newly applied to the severe accident analysis code, SPECTRA, which can conduct dynamic plant evaluation in the different severe accident conditions of nuclear reactors, to develop an evaluation methodology for typical external hazards. In the DSPN-CMMC-SPECTRA coupled frame, the latest safety functions of the plant components/systems can be stochastically determined by the DSPN-CMMC grounded on the current plant states under continuous hazard and the interaction between the multi-state components/systems; then, SPECTRA can evaluate the following plant state determined by the latest safety function of the components/systems. Therefore, the advantage of this newly developed DSPN-CMMC-SPECTRA frame is having the capability to quantitatively and stochastically evaluate the transient accident progressions that potentially lead to the core damage under the continuous external hazard scenario. As for the preliminary exam on the DSPN-CMMC-SPECTRA frame, one of the typical external hazards of continuous volcanic ashfall is selected in this research. In addition, the numerical investigation of alternative accident management' effects has also been carried out and quantitatively confirmed in this research.
Rodriguez, D.; Rossi, F.
Proceedings of INMM & ESARDA Joint Annual Meeting 2023 (Internet), 9 Pages, 2023/05
Under the MEXT subsidy for the improvement of nuclear security related activities, we will present on the current and future progress of the delayed gamma-ray spectroscopic analysis development. We highlight a paper soon to be released and the plan for finalizing the project goal.
久保 光太郎; Jang, S.*; 高田 孝*; 山口 彰*
Journal of Nuclear Science and Technology, 60(4), p.359 - 373, 2023/04
被引用回数：5 パーセンタイル：78.52(Nuclear Science & Technology)確率論的リスク評価(PRA)は、原子力発電所の安全性を向上させるための重要なアプローチである。しかし、この手法では、複合ハザードのモデル化は困難である。地震に起因した溢水シナリオでは、地震による炉心損傷、溢水による炉心損傷、地震と溢水が組み合わさった炉心損傷といった複数の炉心損傷シーケンスを含んでいる。溢水に係るフラジリティは、溢水がタンクなどの水源から区画に伝播するため、時間依存性を有している。そのため、現実的なリスク評価及び定量化を行うためには、動的リスク評価を用いる必要がある。本研究では、地震,溢水,熱水力シミュレーションを連成させ、複数ハザード間の依存関係を明示的に考慮し、地震起因溢水のリスク評価を行った。特に、福島第一原子力発電所事故を踏まえた安全性向上対策に注目し、システムの耐力に関する感度解析と可搬型ポンプを用いた蒸気発生器代替注水の効果を評価した。我々は、シミュレーションに基づく動的PRA手法の複合ハザード起因のリスクの評価への使用を実証した。
Li, C.-Y.; 渡部 晃*; 内堀 昭寛; 岡野 靖
Journal of Nuclear Science and Technology, 23 Pages, 2023/00
被引用回数：1 パーセンタイル：72.91(Nuclear Science & Technology)For all the nuclear reactor systems, quantitative assessment of the accident management (AM) effects against long-term external hazards became one of the essential issues after the lesson learned from the Fukushima Daiichi Nuclear Power Plant accident. However, the influence from the safety systems' stochastic and dynamic shifting between multiple working states, which is related to the interaction with the adjacent components/systems in general, has not been accounted for yet. Therefore, this research aims to develop a dynamic probability risk assessment tool considering repairable multi-component interdependency for investigating the AM influences on the multi-state safety systems under long-term external hazards. Based on the newly proposed methodology in this research via integrating the Petri Net (PN) model with the continuous Markov chain Monte Carlo (CMMC) method, a framework applying PN-CMMC methodology to a severe accident analysis code, SPECTRA, had been originally constructed. Different AM influences on the multi-state decay heat removal systems against long-term volcanic ashfall were also quantitatively confirmed, indicating that halving the repairing time is more influential in suppressing the core damage frequency than doubling the number of adjacent electricity support systems. Therefore, the PN-CMMC-SPECTRA framework can further assess the uncharted dynamic multi-state concerns, leading to a safer AM strategy.
Tuya, D.; 長家 康展
Nuclear Science and Engineering, 15 Pages, 2023/00
被引用回数：0 パーセンタイル：0.01(Nuclear Science & Technology)局所的な分散低減が必要な問題や深層透過問題に対するモンテカルロ中性子輸送計算では、計算効率を上げるため一般的にインポータンス関数の評価が必要である。本研究では、局所的な分散低減に対するインポータンス関数を評価するための新しい再帰的モンテカルロ法(RMC法)を開発した。新しいRMC法を1次元鉄平板上体系と3次元コンクリート鉄体系の2つのサンプル問題に適用した。新しいRMC法によるインポータンス関数を用いて計算された分散低減パラメータを用いて、これらの問題に対する検出器応答を評価した。分散低減法を用いたモンテカルロ計算によって得られた結果は、分散低減法を用いないモンテカルロ計算の結果とよく一致した。さらに、分散低減法を用いたモンテカルロ計算は、FOM値で比較において、数倍から10
倍の効率化を達成することができた。また、検出器位置の光学的な深さが深くなるほど計算効率が向上することもわかった。
Simanullang, I. L.*; 中川 直樹*; Ho, H. Q.; 長住 達; 石塚 悦男; 飯垣 和彦; 藤本 望*
Annals of Nuclear Energy, 177, p.109314_1 - 109314_8, 2022/11
被引用回数：0 パーセンタイル：0.01(Nuclear Science & Technology)Power distribution plays a significant role in preventing the fuel temperature exceeds the safety limit of 1600
C in high-temperature gas-cooled reactors. The experiment to measure the power distribution in the graphite-moderated system was carried out with the Very High Temperature Reactor Critical Assembly facility. In the previous study, the power distribution in the VHTRC was calculated using a nuclear design code system based on diffusion calculation. The results showed a maximum discrepancy of up to 20 between the experiment and calculated values in the axial direction. The large discrepancy occurred near the boundary of fuel and reflector regions. This study describes the evaluation results of pin-wise power distribution of the VHTRC with the Monte Carlo MVP3 code. The calculation results were in good agreement with the measured results. In the axial direction, the discrepancy was less than 1 around the boundary of fuel and reflector regions.
Omer, M.; 静間 俊行*; 羽島 良一*; 小泉 光生
第43回日本核物質管理学会年次大会会議論文集(インターネット), 3 Pages, 2022/11
Gamma-rays originated from laser Compton scattering (LCS) are convenient photon sources for nondestructive interrogation of nuclear materials. LCS can be used with nuclear resonance fluorescence (NRF) and X-ray fluorescence (XRF), the two of which are considered photon-based active interrogation techniques. However, an accurate estimation of the incident LCS 
-ray flux is crucial. The -ray flux is customarily measured using high purity germanium (HPGe) detectors, usually calibrated using standard point-like radioactive -ray sources. These standard sources are entirely different from LCS beams in terms of detection geometry. Therefore, the calibration process must be corrected to meet the LCS beam conditions. Here, we demonstrate how to implement the required corrections and provide experimental validation of these corrections.
Zheng, X.; 玉置 等史; 高原 省五; 杉山 智之; 丸山 結
Proceedings of Probabilistic Safety Assessment and Management (PSAM16) (Internet), 10 Pages, 2022/09
Uncertainty gives rise to the risk. For nuclear power plants, probabilistic risk assessment (PRA) systematically concludes what people know to estimate the uncertainty in the form of, for example, risk triplet. Capable of developing a definite risk profile for decision-making under uncertainty, dynamic PRA widely applies explicit modeling techniques such as simulation to scenario generation as well as the estimation of likelihood/probability and consequences. When quantifying risk, however, epistemic uncertainties exist in both PRA and dynamic PRA, as a result of the lack of knowledge and model simplification. The paper aims to propose a practical approach for the treatment of uncertainty associated with dynamic PRA. The main idea is to perform the uncertainty analysis by using a two-stage nested Monte Carlo method, and to alleviate the computational burden of the nested Monte Carlo simulation, multi-fidelity models are introduced to the dynamic PRA. Multi-fidelity models include a mechanistic severe accident code MELCOR2.2 and machine learning models. A simplified station blackout (SBO) scenario was chosen as an example to show practicability of the proposed approach. As a result, while successfully calculating the probability of large early release, the analysis is also capable to provide uncertainty information in the form probability distributions. The approach can be expected to clarify questions such as how reliable are results of dynamic PRA.
Omer, M.; 静間 俊行*; 羽島 良一*; 小泉 光生
Radiation Physics and Chemistry, 198, p.110241_1 - 110241_7, 2022/09
被引用回数：2 パーセンタイル：53.91(Chemistry, Physical)In beam geometries where a directed -ray beam hits the surface of a coaxial high purity germanium detector (HPGe), the detector efficiency is sensitive to the position where -rays initially hit the detector surface because the structure of the detector is nonuniform. This may cause inaccuracy of the detector efficiency when measured using standard sources that are point-like sources emitting -rays isotropically. Obtaining a precise estimation of the full energy peak efficiency of the coaxial HPGe detector in the beam geometry for on-axis and off-axis measurements requires a Monte Carlo simulation. We performed a Monte Carlo simulation that calculates the detector efficiency in the beam geometry. The effects of the off-axis distance and -ray beam size on the efficiency are quantitatively analyzed. We found that the intrinsic efficiency in the beam geometry is maximized when the beam hits the detector at specific off-axis distances. Our Monte Carlo calculations have been supported by a nuclear resonance fluorescence experiment using laser Compton scattering -ray beam.
久保 光太郎; 藤原 啓太*; 田中 洋一; 白田 勇人*; 荒毛 大輔*; 内山 智曜*; 村松 健*
Proceedings of 29th International Conference on Nuclear Engineering (ICONE 29) (Internet), 8 Pages, 2022/08
福島第一原子力発電所の事故後、外部事象、特に地震や津波に対する確率論的リスク評価(PRA)の重要性が認識された。日本原子力研究開発機構では、地震PRAのための解析手法として、DQFM(direct quantification of fault tree using Monte Carlo simulation)法を開発してきた。DQFMは、地震応答と機器の耐力に関する適切な相関行列が与えられた場合、解析解や多次元数値積分を用いてミニマルカットセット確率を求める方法では困難であるANDゲートやORゲートで接続された機器の相関損傷の影響をフォールトツリーで考慮することが可能であり、その有用性が示されている。このDQFMの計算時間を短縮することは、規制機関や事業者が実施するPRAにおいて、多数の解析が可能になる。そこで、本研究では、準モンテカルロサンプリング,重要度サンプリング及び並列計算の3つのアプローチで予備的検討を行い、計算効率の向上を図った。具体的には、加圧水型原子炉の簡易的なPRAモデルに対する、DQFMによる条件付炉心損傷確率の計算に準モンテカルロサンプリング,重要度サンプリング,並列計算を適用した。その結果、準モンテカルロサンプリングは、仮定した中及び高地震動レベルで有効であり、重要度サンプリングは低地震動レベルで有効という結果が示された。また、並列計算により、実用的な不確実さ解析及び重要度解析が実施可能であることが示された。これらの改良を組合わせてPRAコードに実装することで、計算の大幅な高速化が期待でき、リスク情報を用いた意思決定におけるDQFMの実用化の見通しを得た。
Li, C.-Y.; 渡部 晃*; 内堀 昭寛; 岡野 靖
第26回動力・エネルギー技術シンポジウム講演論文集(インターネット), 4 Pages, 2022/07
シビアアクシデントに至る可能性のある事故シナリオを同定し、その発生頻度を評価することは重要な課題である。本研究ではナトリウム冷却高速炉を対象とし、時間依存や事象の相互依存性を考慮できる動的PRA評価手法の確立を目指す。具体的には過酷事故解析コードSPECTRAに対して新たに連続マルコフ連鎖モンテカルロ(CMMC)を適用し、外部ハザードに対する解析手法を開発する。現在、崩壊熱除去系における空気冷却器のフォルトツリーモデルをCMMCとして実装し、火山降灰に対するプラント過渡特性の試行解析が終了した。
-eigenvalue perturbation in Monte Carlo calculationTuya, D.; 長家 康展
Annals of Nuclear Energy, 169, p.108919_1 - 108919_9, 2022/05
被引用回数：1 パーセンタイル：31.61(Nuclear Science & Technology)核分裂性物質の系に加えられる変化が固有値に与える影響を評価するには、摂動論と言われる特別な方法が必要である。この研究では、相関サンプリング(CS)法と随伴重み法に基づく繰り返し核分裂確率(IFP)を組み合わせることにより、近似を用いない厳密な摂動理論に基づく随伴重み相関サンプリング(AWCS)法を開発した。この厳密な摂動理論に基づいて開発されたAWCS法は、小さな摂動に対しても非常に小さい不確かさを与えるCS法の利点と連続エネルギモンテカルロ法に対して安定した結果を与えるIFP法に基づく随伴重み法の利点を有しており、摂動計算のための新しい正確な方法を与えるものである。開発されたAWCS法の検証のために行ったGodiva及び単純化されたSTACY数密度摂動問題の解析結果は、参照計算結果とよく一致をした。
Riyana, E. S.; 奥村 啓介; 坂本 雅洋; 松村 太伊知; 寺島 顕一
Journal of Nuclear Science and Technology, 59(4), p.424 - 430, 2022/04
被引用回数：1 パーセンタイル：16.35(Nuclear Science & Technology)Modification of the Monte Carlo depletion calculation code OpenMC was performed to enable the depletion calculation of the subcritical neutron multiplying system. With the modified code, it became possible to evaluate the quantity of short half-life fission products from spontaneous and induced fissions in the subcritical system. As a preliminary study, it was applied to the fuel debris storage canister filled with nuclear materials and spontaneous fission nuclides. It was confirmed that the code could successfully provide a quantity of short half-life FPs over time and provide the relationship between the activity ratio of Kr-88 to Xe-135 and effective neutron multiplication factor of the canister.
